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ABSTRACT 

In this video we highlight the application of Computational 
Geometry to our understanding of the formation and dy- 
namics of the Cosmic Web. The emergence of this intricate 
and pervasive weblike structure of the Universe on Mega- 
parseifj] scales can be approximated by a well-known equa- 
tion from fluid mechanics, the Burgers' equation. The so- 
lution to this equation can be obtained from a geometrical 
formalism. We have extended and improved this method by 
invoking weighted Delaunay and Voronoi tessellations. The 
duality between these tessellations finds a remarkable and 
profound reflection in the description of physical systems in 
Eulerian and Lagrangian terms. 

The resulting Adhesion formalism provides deep insight 
into the dynamics and topology of the Cosmic Web. It 
uncovers a direct connection between the conditions in the 
very early Universe and the complex spatial patterns that 
emerged out of these under the influence of gravity. 
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1. THE COSMIC WEB 

The Cosmic Web is the fundamental spatial organization 
of matter on scales of a few up to a hundred Megaparsec. 
Galaxies, intergalactic gas and dark matter exist in a wispy 

1 The main measure of length in astronomy is the parsec. 
Technically a parsec is the distance at which we would see 
the distance Earth-Sun at an angle of 1 arcsec on the sky. 
One Megaparsec is 3.26 million light-years; this is approxi- 
mately the distance between our Milkyway and the nearest 
other galaxy: Andromeda. 
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weblike arrangement of dense compact clusters, elongated 
filaments, and sheetlike walls, amidst large near-empty void 
regions. The filaments are the transport channels along 
which matter and galaxies flow into massive high-density 
cluster located at the nodes of the web. The weblike net- 
work is shaped by the tidal force field accompanying the 
inhomogeneous matter distribution. 

Structure in the Universe has risen out of tiny primordial 
(Gaussian) density and velocity perturbations by means of 
gravitational instability 8, 11 . The large-scale anisotropic 
force field induces anisotropic gravitational collapse, result- 
ing in the emergence of elongated or flattened matter config- 
urations. According to Zeldovich, the collapse of a primor- 
dial cloud (dark) matter passes through successive stages, 
first assuming a flattened sheetlike configuration as it col- 
lapses along its shortest axis. This is followed by a rapid 
evolution towards an elongated filament as the medium axis 
collapses and, if collapse continues along the longest axis, 
may ultimately produce a dense and compact cluster or halo 
[9]. The hierarchical setting of these processes, occurring si- 
multaneously over a wide range of scales complicates the 
picture considerably. 

2. ZELDOVICH & ADHESION MODELS 

The simplest model that describes the emergence of struc- 
ture and complex patterns in the Universe is the Zeldovich 
Approximation (ZA) [11115]. In essence, it describes a ballis- 
tic flow, driven by a constant (gravitational) potential. The 
resulting Eulerian position x(t) at some cosmic epoch t is 
specified by the expression, 

x(t) =q + D(t) u (q), 

where q is the initial "Lagrangian" position of a particle, 
D(t) the time-dependent structure growth factor and uq — 
— V;$o its velocity. The nature of this approximation may 
be appreciated by the corresponding source-free equation of 
motion, 

du , _ . 

— + (it ■ V x )u = 0. 

The use of ZA is ubiquitous in cosmology. One major ap- 
plication is its key role in setting up initial conditions in 
cosmological N-body simulations. Of importance here is its 
nonlinear extension in terms of Adhesion Model. 




in fluid mechanics. The Voronoi diagram represents vol- 
ume. A vertex of the Delaunay triangulation represents 
a mass-element. Its dual Voronoi cell is the volume this 
mass-element occupies. As time moves on, the particle is 
expanding in a struggle for dominance over its neighbours, 
a struggle that is lost once the vertex becomes redundant 
(in the sense that it no longer contributes to the Delaunay 
triangulation). Physically the particle is then trapped in a 
larger structure. In general particles first stream into walls 
(objects of co-dimension 1), moving down into objects of 
increasing co-dimension. [4] [5]. 



Figure 1: Convex hull of the modified potential and 
the resulting weighted Delaunay triangulation. 





Figure 2: Convex conjugate of the modified poten- 
tial and the resulting weighted Voronoi diagram. 



The ZA breaks down as soon as self gravity of the form- 
ing structures becomes important. To 'simulate' the effects 
of self gravity, Gurbatov et al. [2] [3] included an artificial 
viscosity. This results in the Burgers' equation, 

du 2 
— + [u ■ \/ x )u = uV x u, 

a well known PDE from fluid mechanics. This equation has 
an exact analytical solution [6], 



u(x, D) 



Jdq^ exp[G/(»] 



where G(x,q,D) = $o(g) 
v — > 0, the solution is 



fdqexp[G/(2v)] ' 

- (a: - q) 2 /(2D). In the limit of 



§(x, D) = max 

<3 



$o(g) 



if 



2D 



This leads to a geometric interpretation of the Adhesion 
Model. The solution follows from the evaluation of the con- 
vex hull of the velocity potential modified by a quadratic 
term [10]. We found that the solution can also be found 
by computing the weighted Voronoi diagram 7 of a mesh 
weighted with the velocity potential. It is of special im- 
portance that this computation is done on a periodic space. 
This is essential to treat boundaries correctly. 

3. ADHESION & TESSELLATIONS 

The duality between the weighted Voronoi diagram and 
the weighted Delaunay triangulation has a deep connection 
with the concepts of Eulerian and Lagrangian coordinates 



4. THIS VIDEO 

We use the Computational Geometry Algorithms Library 
[T] for our key geometric computations: weighted Delaunay 
triangulations in both 2D and 3D (at this point, the periodic 
triangulation package of CGAL only offers 3D periodic non- 
weighted Delaunay triangulations). 

This video takes us from the initial conditions, through 
both the Zeldovich and Adhesion models, to a display of 
how the topology of the Cosmic Web depends directly on 
the starting point. 
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